Walker et al. BMC Genomics 2014, 15:552 
http://www.bionnedcentral.conn/1471 -21 64/1 5/552 



RESEARCH ARTICLE Open Access 



Genome-wide identification of the Fermentome; 
genes required for successful and timely 
completion of wine-like fermentation by 
Saccharomyces cerevisiae 

Michelle E Walker^'^, Trung D Nguyen^'^, Tommaso Liccioli^'^, Frank Schmid^'^, Nicholas Kalatzis^'^ 
Joanna F Sundstrom^'^, Jennifer M Gardner^'^ and Vladimir Jiranek^'^"" 



Abstract 

Background: Wine fermentation is a harsh ecological niche to which wine yeast are well adapted. The initial high 
osmotic pressure and acidity of grape juice is followed by nutrient depletion and increasing concentrations of 
ethanol as the fermentation progresses. Yeast's adaptation to these and many other environmental stresses, enables 
successful completion of high-sugar fermentations. Earlier transcriptomic and growth studies have tentatively 
identified genes important for high-sugar fermentation. Whilst useful, such studies did not consider extended growth 
(>5 days) in a temporally dynamic multi-stressor environment such as that found in many industrial fermentation 
processes. Here, we identify genes whose deletion has minimal or no effect on growth, but results in failure to achieve 
timely completion of the fermentation of a chemically defined grape juice with 200 g total sugar. 

Results: Micro- and laboratory-scale experimental fermentations were conducted to identify 72 clones from -5,100 
homozygous diploid single-gene yeast deletants, which exhibited protracted fermentation in a high-sugar medium. 
Another 21 clones (related by gene function, but initially eliminated from the screen because of possible growth 
defects) were also included. Clustering and numerical enrichment of genes annotated to specific Gene Ontology (GO) 
terms highlighted the vacuole's role in ion homeostasis and pH regulation, through vacuole acidification. 

Conclusion: We have identified 93 genes whose deletion resulted in the duration of fermentation being at least 20% 
longer than the wild type. An extreme phenotype, 'stuck' fermentation, was also observed when D0A4, NPTl, PLCl, 
PJK2, SINS, SSQl, TPSl, TPS2 or ZAPl were deleted. These 93 Fermentation Essential Genes (FEG) are required to 
complete an extended high-sugar (wine-like) fermentation. Their importance is highlighted in our Fermentation 
Relevant Yeast Genes (FRYG) database, generated from literature and the fermentation-relevant phenotypic 
characteristics of null mutants described in the Socchoromyces Genome Database. The 93-gene set is collectively 
referred to as the 'Fermentome'. The fact that 10 genes highlighted in this study have not previously been linked 
to fermentation-related stresses, supports our experimental rationale. These findings, together with investigations 
of the genetic diversity of industrial strains, are crucial for understanding the mechanisms behind yeast's response 
and adaptation to stresses imposed during high-sugar fermentations. 
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Background 

Saccharomyces cerevisiae strains are widely used for pro- 
duction of alcoholic beverages such as wine, beer, sake, 
as well as bioethanoL The selection of yeast strains for 
efficient fermentation performance in these industrial 
processes has typically focused on attributes such as pre- 
dictable fermentation at the relevant process tempera- 
tures, desired fermentation vigour and extent of sugar 
attenuation with efficient conversion to ethanol. Attri- 
butes deemed important for wine fermentation include 
retention and enhancement of varietal (grape derived) 
fruit characters and production of desirable flavour and 
aroma compounds [1]. The numerous subtle differences 
in fermentation traits between the various industrial 
yeast strains are reflected in the overall genetic diversity 
found in this species [2]. Nevertheless, industrial yeast 
share the ability to grow and ferment in high sugar me- 
dia, sensing, responding and adapting to the extreme 
and changing conditions. Such conditions are imposed 
by anaerobic conditions, high concentrations of sugar 
and organic acids (resulting in high osmotic pressure), 
toxins and inhibitors, at times low pH and assimilable 
nitrogen levels, extended fermentation times and in- 
creasing concentration of ethanol or other inhibitors. 
Our understanding of the cellular mechanisms behind 
yeasts adaptation to the temporal exposure to multiple 
environmental stresses in such fermentations is limited 
despite extensive studies in the past decade. 

Extensive phenomic' studies have been undertaken 
with laboratory yeast collections, comprised of individ- 
ual known single gene deletion mutants (deletants), 
whereby the phenotype of such deletants is analysed to 
determine the genes associated with tolerance to a spe- 
cific condition. Several studies have looked at tolerance 
or growth sensitivity to singular conditions related to 
high sugar fermentation, including high osmolarity (spe- 
cifically sucrose [3] and glucose [4]), anaerobic growth 
[5], oxidative stress [6], tolerance to high pressure and 
low temperatures [7], ethanol [8-10] and acetic acid 
[11]. These studies typically did not require the comple- 
tion of fermentation and were conducted under aerobic 
conditions usually for 1-3 days only. 

This study is unique in that it sought to identify genes 
required by yeast to not only grow in but to complete 
fermentation of a high sugar medium, that is a che- 
mically defined grape juice, wherein multiple stresses 
would be experienced simultaneously or sequentially. 
The non-availability of a deletion library in a protot- 
rophic wine yeast background necessitated the use of 
the laboratory yeast library, derived from S288c [12]. 
Although these lab yeast strains are reported to cope 
poorly with high osmotic pressure as found in grape 
juice [13], we and others have shown that the S288c de- 
rived diploid BY4743 and haploids BY4741 (data not 



shown) and BY4742 [13], are able to complete fermenta- 
tion in high sugar media [14,15], provided their auxo- 
trophic requirements are met. 

Previous studies have identified genes essential for 
growth under particular conditions relevant to fermenta- 
tion, but not the entire fermentation process itself. Suffi- 
cient biomass is critical for fermentation performance; if 
there is too little, fermentation becomes severely pro- 
tracted or may even fail to complete. Adequate biomass 
does not however guarantee fermentation completion, es- 
pecially in high sugar media. Furthermore, whilst transcrip- 
tomic experiments highlighted genes whose expression 
was significantly altered during fermentation, these genes 
are not necessarily required for fermentation reliability. 

This report relates to the effect of gene deletions on 
yeast s ability to complete extended high sugar fermenta- 
tion. Those deletants whose maximum biomass was at 
least 70% of the parent but which exhibited protracted 
(or stuck) fermentation were taken to highlight genes 
fundamental for successful fermentation in high sugar 
conditions, such as in wine production. The clustering 
of these into specific groupings i.e. Gene Ontology (GO) 
terms provides insight into the many cellular mecha- 
nisms behind yeast s adaptation to fermentation stresses. 

Results and discussion 

A micro-fermentation screen to identify genes required 
for successful fermentation in a high sugar medium 

In this study (outlined in Figure 1), a high throughput 
fermentation screen was developed to identify those ge- 
nes required for high sugar (i.e. akin to wine) fermenta- 
tion, whereby deletion of the corresponding genes would 
result in protracted or arrested fermentation. Whereas 
previous screens have focused on growth assays using 
standard laboratory media such as YPD supplemented 
with a stress agent (e.g. ethanol, sugar, oxidizing com- 
pounds), our approach was to test the ability of yeast to 
complete extended fermentation, albeit on a micro-scale. 
The ability to undergo a wine fermentation requires the 
yeast to respond and adapt to a dynamic and temporal 
exposure to multiple stresses in the 'juice' environment. 
Specifically, upon inoculation the cell population must 
be able to grow to a high density whilst exposed to high 
osmotic conditions (imposed by the sugars and organic 
acids in the juice) and the rapid depletion of oxygen and 
nitrogenous nutrients, sterols and vitamins [16]. With 
the bulk of fermentation being conducted by cells in sta- 
tionary phase, long-term adaptation is paramount to the 
cells tolerance to increasing levels of ethanol, as a con- 
sequence of sugar uptake and catabolism. 

Chemically defined grape juice medium (CDGJM_200) 
containing 200 g L"^ sugar (as equimolar amounts of 
glucose and fructose) was chosen for the high sugar 
medium. Whilst fructose and glucose are found in grape 
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Figure 1 Outline of fermentation study. Schematic of the fermentation screen (micro-fermentations; 0.6 mL and laboratory-scale; 100 mL 
fermentations) and number of gene deletion mutants identified at each stage as having protracted fermentation. Final evaluation using 100 mL 
fermentations identified 72 candidate genes. In silico analysis of fermentation relevant data (SGD_FRYG and LIT_FRYG databases) and GO analysis 
of genes functionally related to 72 candidate genes. 21 additional gene deletants were identified as resulting in protracted fermentation. 
Together, the 93 genes are referred to as Fermentation Essential Genes or the Termentome'. FD denotes fermentation duration. 



must [17] and 200 g L"^ of sugars is considerably higher 
than the 20 g L"^ used in laboratory media, this con- 
centration is at the lower end of typical concentrations 
found in grape must (240-300 g L"^ [18]). This com- 
promise was made due to the poorer fermentative ability 
of laboratory strains in comparison to wine yeast strains 
[19], in particular the reported sensitivity of S288c deriv- 
atives under high osmotic pressure [13]. The assimilable 
nitrogen content (487 mg N L"^) was sufficient for com- 
plete fermentation of 200 g L"^ of sugar and greater than 
the minimum requirement by yeast (140 mg N L"^ [15]) 
and also addressed the auxotrophic requirements of the 
diploid strain, BY4743, in which the yeast homozygous 
gene deletion collection was generated [12]. At the study's 



commencement the lack of availability of a prototrophic 
yeast deletion collection in either a laboratory yeast [20] 
or wine yeast genetic background [21] meant that fermen- 
tation without the addition of amino acids matched to the 
auxotrophic requirements of BY4743 could not be investi- 
gated. Such an investigation is nonetheless of interest to 
determine the genetic basis for differences in the nitrogen 
efficiency of yeast [15]. 

The initial screen was performed in duplicate 0.6 mL 
fermentations at 28°C over 7 days in uracil-supplemented 
CDGJM_200. This medium supported growth of the par- 
ental strain BY4743, with all glucose and fructose being 
consumed within 144 h (data not shown). The individual 
fermentation performance of each of the -5100 yeast 
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clones of the homozygous diploid deletion collection 
was compared with BY4743 (data not shown; Figure 1). 
The extent of fermentation was defined by enzymatic de- 
termination of the residual glucose and fructose content 
after 180 h and was considered complete when the total 
residual sugar was <2.5 g L"^. The additional 36 h, repre- 
senting an extra 25% in fermentation duration, was chosen 
to allow for growth discrepancies due to the potential for 
minor differences in inoculation rates. 336 gene deletion 
mutants were identified as having failed to complete fer- 
mentation after 180 h (Additional file 1). 

In order to eliminate the possibility that incomplete fer- 
mentation by these 336 deletants was due to impaired 
growth, micro-fermentations were repeated with the add- 
itional estimation of growth (OD 600 nm) over the dur- 
ation of the fermentation (Figure 1 and Additional file 1). 
Residual sugars in the subsequent experiment were ana- 
lysed at 137 h (BY4743 residual sugar = 2.48 g L"^), and 
only deletants with optical densities (ODeoo) comparable 
to BY4743 (i.e. > 70%) were chosen for further investiga- 
tion. In such candidates, the corresponding gene deletion 
was considered to have no or little effect on growth and to 
more specifically affect fermentation. As such, 101 dele- 
tants were chosen as potential candidates dysfunctional in 
fermentation under high sugar conditions. 



Evaluation of deletants in laboratory scale (100 mL) 
fermentation 

To test that the fermentation phenotype of the 101 de- 
letants was reproducible in a more controlled environ- 
ment, the performance of these was compared with the 
parent strain in larger, 100 mL fermentations (Figure 1). 
Fermentations conducted at this scale allow greater con- 
trol over parameters which influence experimental repro- 
ducibility and fermentation outcome, such as inoculum 
preparation, maintenance of anaerobiosis and the inclu- 
sion of biological triplicates. The relative fermentation 
duration (RFD) of deletant versus parent was recorded in 
the experimental data (Additional file 1). Seventy two 
deletants had protracted fermentation, taking 20% or lon- 
ger to complete fermentation than the parent. Eight of 
these were not only protracted but did not complete fer- 
mentation (Figure 2) and the fermentations were deemed 
stuck'. 

The remaining 29 mutants exhibited fermentation du- 
rations similar to the parent (14 deletants, RFD = 1.0; 13 
deletants, RFD = 1.1; 2 deletants, RFD <1). These findings 
indicate that although there was a 29% over-estimation of 
genes affecting fermentation in the second screen (where 
optical density was measured), 71% of the clones were stiU 
identified as affecting fermentation. 
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Figure 2 Comparative fermentations with parent yeast BY4743 and nine stucic mutants in CDGJIVI_200. Fermentations (100 mL) were 
performed in triplicate in CDGJM_200 at 28°C wliereby tine deletants were compared directly with the parental strain BY4743. Growth was 
monitored as optical density at 600 nm (open symbols). Sugar consumption was monitored enzymatically and reported as total or individual 
sugars (solid symbols). The data are arranged in three rows (A-C) of three plots each (left to right). The first plot in each row depicts total residual 
sugar. Glucose-only and fructose-only information is shown to the right of the corresponding total sugar plot. The nine deletion strains (deletants) 
which result in arrested or 'stuck' fermentation are shown in rows A (Anptl, Aplcl, Assql), B (Adoa4, AsinS, Azapl) and C (Aptk2, Atpsl, Atps2). 
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The 72 deletants identified as being protracted in fer- 
mentation in a high sugar medium (CDGJM_200) were 
screened as part of the original 336 (Figure 1) in a 
micro-fermentation of low sugar medium (CDGJM_20) 
containing 20 g L"^ (Additional files 1 and 2). The 72 
deletants could be separated into two distinct groups 
based on two time points at which residual sugar was 
measured during fermentation: 32 h (fermentation com- 
pleted by the parent) and 56 h (175% duration). At 32 h, 
the cut-off of 2.5 g L"^ residual sugar was used; 43 dele- 
tants were found to have completed fermentation, with 
29 deletants having >2.5 g L"^ sugar. At 56 hours all 
the deletants except Atkll had completed fermentation. 
TKLl was retained in the 72 gene dataset based on the 
delay in fermentation when examined in 100 mL scale 
using CDGJM_200. 

Use of the FRYG databases to identify genes which 
modulate fermentation - the 'Fermentome' 

Two FRYG databases (Fermentation T^elevant Yeast Genes) 
were separately compiled from previously published 
fermentation relevant studies and the Saccharomyces 
Genome Database (Additional file 3). The databases were 
complementary; the first (LIT_FRYG) was derived from 
relevant literature, examining single stress conditions 
related to fermentation. The second (SGD_FRYG), was 
compiled from fermentation relevant phenotype terms 
specifically selected from the more general SGD pheno- 
type terms database, whereby only null mutants (dele- 
tants) with a specific phenotypic response, i.e. decreased 
resistance or increased sensitivity were selected. The 72 
gene dataset (representing the remaining candidates fol- 
lowing 100 mL fermentations) was compared to the two 
FRYG databases in order to determine which were previ- 
ously reported (Figure 1). The SGD_FRYG database was 
more inclusive with 65 out of the 72 genes being previ- 
ously reported in related conditions whilst only 55 out of 
the 72 genes were identified in the LIT_FRYG database. 

GO analysis of the identified 72 fermentation essential 
genes highlighted specific biological processes and more 
importantly, the other gene members for each GO term, 
that are likely to be important for fermentation. To en- 
sure that these related genes were not incorrectly ex- 
cluded as 'false negatives' they were re-examined here in 
a follow up screen. Of the 29 mutants examined, 21 gave 
protracted fermentation in CDGJM_200. 

One of these, (Atpsl) did not finish fermentation and 
was considered stuck' (Figure 2 and Additional file 2). 
The 21 genes were cross-referenced to both FRYG data- 
sets (Additional file 3). The genes absent from the data- 
bases are shown in Table 1. 

In summary, 93 genes were identified as being required 
for fermentation, whereby deletion of the gene resul- 
ted in the yeast undergoing protracted fermentation 



Table 1 Fermentation essential genes shown to be absent 
in either LIT FRYG or SGD FRYG databases 



FEGs absent in 
LIT FRYG 



FEGs absent 
in SGD_FRYG 



FEGs absent in LIT_FRYG 
and SGD_FRYG 



DUFl {YOL087Q 
GPA2 {YER020W} 
NHXl {YDR456W) 
OPIl {YHL020Q 
PUGl {YER185W} 
RBL2 {YOR265W) 
RTT103 {YDR289Q 
TOM! {YDR457W} 
VAC8{YEL013W} 
YLL007C 
ZAPl {YJL056Q 
GPRl {YDL035Q 
OPIl {YHL020Q 



HRKl {YOR267Q 
CYK3 {YDLn7W} 



CCZl {YBR131W} 
ATG7 {YHR171W} 
HXKl {YFR053Q 
PUGl {YER185W} 
RXT3 {YDL076Q 
YFL012W 
CISl {YDR022Q 
SNX4 {YJL036W) 



The 93 FEG dataset identified after the fermentation screen and GO analysis 
depicted in Figure 1, was compared with the LIT_FRYG and SGD_FRYG 
databases (see Additional file 3 to determine which genes were absent and 
thus unique to this study). 

in CDGJM_200 (Additional file 2). Together, the 93 genes, 
are referred to as the FEG dataset (for Fermentation Es- 
sential Genes) and comprise the laboratory yeast Termen- 
tome'; genes that are essential for the timely completion of 
fermentation. The Termentome' is taken to represent the 
genes/processes required for yeast to sense and respond 
to the multiple stresses of the juice environment. These 
cellular processes enable yeast to grow in grape juice and 
complete fermentation over the extended period of 4-12 
days typical of a wine fermentation. The experimental de- 
sign however did not allow for the identification of genes 
whose deletion shortened fermentation, as no time point 
was taken prior to the parent finishing. Independent stud- 
ies have identified genes which affect carbon and nitrogen 
flux and thus positively modulate fermentation, although 
the mechanisms remain unclear. For example the genes 
NGRl and GID7, upon disruption, resulted in enhanced 
catabolism of sugar in a wine and laboratory strain during 
growth in CDGJM_200 with limiting nitrogen [22]. Also, 
deletion of PDAl, encoding the pyruvate dehydrogenase 
Ea subunit, was shown to alter carbon flux during anaer- 
obic fermentation, resulting in shortened fermentation, 
after an extended lag phase (12 h), which the authors 
suggested was an acclimatisation to the fermentation 
conditions [23]. 

The studies chosen to build the LIT_FRYG dataset in- 
cluded those examining growth in response to a single 
stress relevant to wine fermentation [3-11,24-26] as well 
as a transcriptomic study of the wine strain Vin 13 during 
fermentation of a Riesling juice [27]. Our dataset repre- 
sents only a small percentage of the genes identified in the 
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above studies, as shown in LIT_FRYG (Additional file 3), 
relating to growth sensitivity as a measure of stress toler- 
ance. Not surprisingly, the largest representations relate to 
hypertonic conditions (21.4%) [24], high glucose (17%) [4], 
ethanol (8.3% [9], 18.4% [9], 14.9% [10]), anaerobic 
growth (13.6% [26], 8.3% [5]) and sucrose (9.9%) [3] - see 
Additional file 3. The importance of experimental de- 
sign is highlighted by the fact that our study con- 
tained only 1.79% of genes described as Fermentation 
Stress Response (FSR) genes [27]. Although Marks and 
co-workers [27] looked at gene expression in response to 
increasing ethanol concentration during fermentation, 
small changes in transcript levels of FSR may not eventu- 
ate in altered fermentation kinetics or outcome - a more 
likely observation with complete deletion as per our study. 
A similar conclusion was drawn from data relating to 
growth rate and stress tolerance [28], whereby the au- 
thors examined the yeast deletion collection in con- 
tinuous culture. 

Eighteen of the 93 genes, which resulted in protracted 
fermentation when deleted, were annotated to vacuole 
function, specifically the vacuolar ATPase complex. 
Within this dataset, 6 genes were identified by relation- 
ship rather than in the original fermentation screen (72 
gene dataset) due to clonal errors within the library used 
in the study. New clones were sourced (Invitrogen) and 
verified by sequencing prior to analysis (100 mL scale). 
The occurrence of 10 of these genes in the freeze thaw 
stress dataset {VMAl, VMA2, VMA3, VMA4, VMA7, 
VMA8, VMAll, VMA13, VPHl and VPH2) [29] indicates 
that cellular mechanisms required for tolerance to 
freeze-thaw stress [29] may also be needed for successful 
fermentation. By comparison of our results to previous 
studies (see Additional file 3), we propose that the role 
of the vacuole, specifically maintenance of vacuolar acid- 
ification (discussed later) is a key mechanism of cellular 
response and adaptation during the transitory and dy- 
namic stress conditions imposed by acidic pH, high os- 
molarity (sugar, organic acids), increasing ethanol content 
and elevated temperatures during fermentation. 

The association between cross-tolerance to multiple 
stresses and fermentation outcome, as well as fermenta- 
tive growth and stress tolerance was examined. Genes 
common to the FEG dataset and selected datasets from 
the SGD_FRYG database can be highlighted (Figure 3A-F). 
Here only 83 fermentation essential genes are com- 
pared as the additional 10 (of the 93) unique to our 
study are necessarily excluded. Of particular interest 
was the large number of genes, when deleted, were 
associated with decreased resistance to multiple stresses: 
ethanol, acidic and hyper-osmotic conditions and heat 
(for example D0A4 {YDR069C), PLCl {YPL268W), 
TPS2 {YDR074W), and SSQl {YLR369W) as shown in 
Additional file 3). 




Figure 3 Comparison of fermentation essential genes with the 
SGD_FRYG database: phenotypic response to specific stress 
conditions. Ninety three fermentation essential genes (FEG) were 
identified in tliis study. Ten of tlie genes liave not been previously 
associated with fermentation related phenotypes and are unique to 
this fermentation study. The Venn diagrams therefore depict the 
number of genes within the 83 FEG dataset which upon deletion, 
leads to either increased sensitivity or decreased resistance to 
conditions relevant to fermentation, such as anaerobic growth, 
ethanol toxicity, osmotic and oxidative stress and temperature 
extremes. 

I J 

Yeast deletants encoding the vacuolar H"^ ATPase 
(vma mutants) were often associated with cross tolerance 
towards several stress conditions but not necessarily 
fermentative growth. Surprisingly, none of the 83 gene de- 
letions were associated with decreased utilisation of nitro- 
gen and fermentative metabolism (Additional file 3). 

Classification of genes Into specific gene ontology 
categories 

Using GO toolbox [30] and GO module [31] the Fer- 
mentation Essential genes (FEG) were classified into dis- 
tinct groups according to GO terms, based on their role 
in a biological process, function within the cell and cel- 
lular localisation (Additional file 4). The frequency of oc- 
currence within the genome of genes annotated to these 
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GO terms was compared to the frequency of occurrence 
of genes within the FEG dataset annotated to these same 
terms. An increased frequency of occurrence (an enrich- 
ment) of genes in the FEG dataset annotated to a spe- 
cific GO term impUes the importance of that GO cluster 
to fermentation completion. 

A subset of the most enriched GO terms {p = 10"^^ to 
10"^^) includes gene clusters related to vacuolar function 
and cellular pH and ion homeostasis (Figure 4). The 
vacuole has a multifunctional role; trafficking of mem- 
brane proteins to the plasma membrane, proteolytic deg- 
radation and recycling of proteins sent to the vacuole 
[32], sequestration of toxic metal ions, ion homeostasis 
(calcium and phosphate), osmo-regulation (Na"^/lC ac- 
cumulation) and storage of amino acids [33]. Disruption 
of these functions can therefore have far reaching conse- 
quences for a cell grown at high osmolarity and low pH 
conditions. 



Incomplete 'stuck' fermentation as a result of deletion of 
specific genes 

Nine deletants within the 93 gene FEG dataset gave rise 
to incomplete (stuck) fermentations (Figure 2). These 
genes are central to ion homeostasis {PTK2 {YJR059W)), 
SSQl {YLR369W)), NAD recycling (NPTl {YOR209C)), 
signalling {PLCl {YPL268W)), trehalose synthesis {TPSl 
{YBR126CI TPS2 {YDR074W)), transcription {SIN3 
{YOL004W), ZAPl {YGR28SC)), and ubiquitin recycling 
{D0A4 {YDR069C)), With the exception of Atpsl and 
Aptk2, which showed normal growth, the other 7 stuck' 
mutants exhibited reduced biomass (as determined by 
maximum ODeoo)- Our findings are only in part compar- 
able to the SGD database where Atpsl, Atps2, Anptl, 
Aplcl, and Assql were shown to have reduced fermenta- 
tive growth (Additional file 3 and Table 2). Underlying 
these differences is likely the experimental approach 
used. It is evident that final biomass (maximal OD^oo) 




o 
c 
o 
E 

Figure 4 Enrichment of FEG annotated to specific GO terms from SGD using GO ToolBox GO-Stats. The frequency of gene occurrence 
(expressed as a percentage) for individual GO terms was compared between tine 93 FEG dataset and tine entire yeast genome (7168 genes). 
Tliese GO terms were not mutually exclusive. The statistical probability was calculated to determine the enrichment of the specific GO terms, 
expressed as p-values. Numbers in brackets represent the number of genes in FEG annotated to a specific GO term. Data from the yeast genome 
is denoted with □ and the FEG dataset (93 genes) with ■. 
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Table 2 SGD_FRYG fermentation phenotypes associated with nine 'stuck fermentation' mutants 


Acidic pH 


Hyperosmotic 


Oxidative stress 


Cold 


Heat Starvation 


Ethanol 


resistance 


stress resistance 


resistance 


sensitivity 


sensitivity resistance 


resistance 




550/ 


550/ 


550/ 


550/ 


550/ 


TPS2 


TPS2 


TPS2 




TPS2 


TPS2 


D0A4 


D0A4 


D0A4 




D0A4 


D0A4 




PLC! 


PLC! 




PLCl 


PLCl 






PTK2 




PTK2 


PTK2 






NPTl 




NPTl 




TPSl 










TPSl 










SINS 


SINS 


ZAP! 




ZAPl 









Phenotype responses which result in decreased resistance are denoted with a downward (|) arrow, whilst those which result in increased sensitivity are denoted 
with an upward (f) arrow. 



does not necessarily determine whether yeast are able to 
catabolise sugar albeit in an extended time frame (pro- 
tracted) or fail to catabolise all of the sugar (defined as a 
stuck fermentation). Four vma mutants (defective in vacu- 
olar ATPase) that had a similar growth deficiency (plus an 
extended lag phase) were able to catabolise sugar albeit 
more slowly than the parent BY4743 (Figure 5). 

All 9 genes associated with stuck fermentations were 
shown to have increased sensitivity towards multiple 
stresses relevant to fermentation (Table 2 and Additional 
file 3), with D0A4, SSQl and TPS2 being associated with 
five stresses. Interestingly, only SINS, SSQl, TPSl and 
TPS2 genes were associated with response to stress' 
GO:0006950 (Additional file 4) although there was no 
over-representation within the FEG dataset for genes an- 
notated to this GO term. Whilst Aplcl results in de- 
creased glucose utilisation, none of the 9 stuck mutants 
are listed as decreased fermentation metabolism, de- 
creased nitrogen utilisation or decreased anaerobic growth 
(Additional file 3). 

Role of TPSl and TPS2 as protectants against multiple 
stresses 

TPSl (trehalose-6-phosphate synthase) and TPS2 (tre- 
halose-6-phosphate phosphatase) are required for trehal- 
ose biosynthesis and not surprisingly were both identified 
in this study as FEGs. Deletion of either TPSl or TPS2 
leads to complete cessation of sugar utilisation in the 
second half of fermentation i.e. a stuck' fermentation 
(Figure 2). Trehalose, a non-reducing disaccharide, or- 
dinarily accumulates during stationary phase and is 
classified as a storage carbohydrate [34]. It is associ- 
ated with increased survival and cell protection from 
stresses such as heat, ethanol and freezing [35]. Ac- 
cordingly, trehalose protects against lipid peroxidation 
and protein carbonylation resulting from ethanol in- 
duced oxidative stress, which would otherwise impact 
on membrane dynamics and glycolysis, respectively [36]. 



Past studies have shown sensitivity of Atpsl to ethanol 
[10] and oxidative stress [6], and Atps2 to ethanol, heat, 
and NaCl [24], high acidity [25] and oxidative stress 
[6,24]. Conversely, over-expression of these genes results 
in enhanced thermotolerance and improved ethanol re- 
sistance [35-37]. Given the recent link between mitochon- 
drial mutation and ageing [38], and trehalose [36], and 
that fermentation is conducted by metabolically active but 
stationary phase cells, it is likely that the lack of trehalose 
in Atpsl and Atps2 [36] results in a vastly decreased abil- 
ity of the cell to protect against the toxic effect of ethanol, 
leading to increased mitochondrial dysfunction and pre- 
mature ageing and reduced viability during fermentation. 

Role of SIN3 in thermal stress 

Adaptation to elevated temperatures, especially above 
30°C and in the presence of ethanol [39], is fundamental 
to survival and the ability to complete fermentation. The 
primary mechanism is the heat shock response operating 
through chromatin modulation by the Rpd3L deacetylase 
complex, which is recruited to target promoters upon 
heat stress affecting transcription and general metabolic 
processes (reviewed in [40]). SINS is associated with the 
heat shock response; deletion of which resulted in the 
complete arrest of fermentation (Figure 2). Sin3p is a 
component of the Rpd3L complex, along with Sap30p 
and Rxt3p. The latter two genes, when deleted, resulted 
in only protracted fermentation. The Rdp3L complex is 
proposed to facilitate the ability of yeast to tolerate mul- 
tiple stresses including thermal and osmotolerance but 
not oxidative stress [28]. 

SSQ1 and ZAP1 in ion homeostasis 

Osmotic stress, whether salt or sugar induced, is well 
known to have a direct effect on cellular ion homeostasis 
[41]. In response to the external environment, various 
transporters transport cations into the organelles and 
across the plasma membrane, such that cell turgidity is 
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Figure 5 Comparative fermentations with parent yeast BY4743 
and four Avma mutants and Aopil in CDGJIVI_200. 

Fermentations (100 mL) were performed in triplicate in CDGJM_200 
at 28°C wliereby tine deletants were compared directly with the 
parental strain BY4743. Fermentation progress was monitored by 
determining the total residual sugar (glucose and fructose) by 
enzymatic analysis. Growth was monitored as optical density 
(600 nm) 20 h post-inoculation and at intervals throughout 
fermentation. 



maintained [42]. In this study 26 genes (28.9% of FEG 
genes) involved in ion homeostasis were identified as 
having a role in fermentation. 

Two FEGs in particular, SSQl and ZAPl, involved in 
iron and zinc homeostasis respectively, resulted in stuck' 
fermentation (Figure 2). SSQl encodes a mitochondrial 
hsp70-type molecular chaperone that is required for the 
assembly of Fe-S clusters into Fe-S proteins, such as the 
mitochondrial proteins aconitase, Yfhlp and several 
DNA repair enzymes [43]. Iron homeostasis is proposed 
to be highly regulated, with the Fe-S clusters thought to 
function as regulatory sensors to oxidative stress and 
intracellular iron (reviewed in [43]). Yeast strains lacking 
the Fe-S cluster scaffold protein Isulp or chaperone pro- 
tein Ssqlp increase iron sequestration to the mitochon- 
dria, and show increased mitochondrial oxidative damage 



[43]. Transcription of ISUl was altered during fermenta- 
tion in response to increasing external ethanol content 
[27], however, in this study Aisul had comparable fer- 
mentation duration to the parent (microferment 1; see 
Additional file 1). Given the speculation on Fe-S proteins 
acting as iron sensors, decreased formation of these pro- 
teins as observed in Assql [44] may lead to a signalling of 
iron deficiency within the mitochondria. This may cause 
an increase in mitochondrial uptake and decrease in ex- 
port of iron, and consequently mitochondrial iron over- 
load and subsequent oxidative damage. Assql has been 
previously shown to be sensitive to oxidation [44], which 
may explain the fermentation phenotype of Assql in 
our study, in terms of oxidative stress associated with 
fermentation in CDGJM_200 and the low iron con- 
tent (14 [ig L"^) [17,45], which is in the lower range 
(20-330 (ig L"^) of typical Australian Chardonnay 
juices [46]. Given that mitochondrial dysfunction has 
recently been linked with vacuolar acidification in af- 
fecting chronological ageing [38] it is of interest that 
the Assql deletant is not defective in vacuolar acidifi- 
cation ([47], data not shown). 

It appears that zinc homeostasis is also critical to fer- 
mentation, through the transcriptional control of some 
80 genes exerted by the zinc-sensing transcription factor 
Zaplp, (reviewed in [48]). These include genes involved 
in phospholipid biosynthesis, zinc uptake and vacuolar 
storage (detoxification), sulfur metabolism, cell wall func- 
tion and response to oxidation (ROS) [48]. Zaplp has a 
protective role with regards to oxidative stress and dam- 
age. The down-regulation of sulfur metabolism by Zaplp 
in zinc deficient cells is associated with increased pro- 
tection through the redirection of NADPH from sulfur 
metabolism to regenerating reduced peroxiredoxin and 
glutathione associated with antioxidant defence mecha- 
nisms (cited in [48]). The extended and incomplete fermen- 
tation profiles of Azapl (Figure 2) indicate that Zaplps 
role as a core regulatory protein particularly in oxidative 
protection and perhaps vacuolar storage is vital for fermen- 
tation maintenance. 

Influence of intracellular pH on fermentation 

GO analysis of the FEG dataset (Additional file 4), not 
only highlighted the enrichment of 18 genes associated 
with the vacuolar ATPase complex (targeted by 
Vphlp isomer of Vq subunit a) [49-51] but also DBF2 
{YGR092W) [52] and NHXl {YDR4S6W) [53], known to 
be involved with regulation of intracellular pH. The tight 
regulation of intracellular pH, whilst influenced by exter- 
nal pH and nutrient availability, is reflective of most 
physiological processes being pH dependent [48]. Crit- 
ical to this, is the cytosols buffering capacity, originating 
from the acid-base action of metabolites such as ammo- 
nium, phosphate and carbon dioxide, as well as proteins 
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[54]. As such, pH sensing' is associated with many pro- 
cesses including carbon flux through the glycolytic cy- 
cle, coupled reduction-oxidation reactions involved in 
NAD(P)^ cycling as well as phosphatidylinositol phos- 
phates acting as signalling molecules to effect trans- 
cription [54]. 

The main pH regulator in yeast is a P2 type H^ATPase; 
Pmalp, responsible for the active translocation of protons 
across the plasma membrane into the external medium, 
thus maintaining cytosolic pH at neutral [48,55]. The sec- 
ond, is the vacuolar (V) H"^ATPase, a complex of two do- 
mains, Vi (8 subunits) and Vq (6 subunits), which is 
responsible for the active transport of protons from the 
cytosol into the vacuolar lumen [49]. PMAl {YGL008C) is 
an essential gene and so not present in the homozygous 
diploid deletion library used in this study. However, the 
heterozygous diploid {PMAl/Apmal) results in haplo- 
insufficiency, i.e. reduced growth rate in grape juice 
conditions (202 mg L"^ YAN, 130 g L"^ water soluble 
carbohydrate (WSC)) and surprisingly also the reverse, 
haploproficiency (increased growth rate) in nitrogen and 
carbon limiting conditions (118.7 mg L"^ YAN, 2.5 g L"^ 
WSC) [56]. 

Pmalp activation is dependent on a complex inter- 
action between glucose sensors such as SnfSp and Rgt2p, 
and the Gprlp/Gpa2p receptor/G protein-coupled (GPCR) 
complex [57], protein kinase C (Pkclp) [57], phospholipase 
C (Plclp) [58], the proteosome, and interaction with tu- 
bulin (as reviewed in [54]). Several studies conducted 
on glucose starved cells, using specific inhibitors of Plcl 
(3-nitrocourmarin) [59,60] and Pmalp (carbonyl cyanide 
m-chlorophenylhydrazone; CCCP) [61] in conjunction 
with glucose additions (up to 100 mM or 18 g L"^) have 
demonstrated the activation of Pmalp to be mediated by 
Plclp and the GPCR complex. Plclp is proposed to trans- 
duce the glucose signal, through calcium signalling, by 
the induction of calcium uptake into the cell and vacuolar 
transport into the cytoplasm [59,62-64] Gprlp aids 
rapid adaption to extracellular glucose (and not fruc- 
tose), via the cAMP-PKA pathway [65]. The latter ac- 
tivates cAMP signalling via a glucose phosphorylation 
mechanism [66]. 

In this study genes associated with hexose sensing were 
identified; GPA2 {YER020W) and GPRl {YDL03SC)), to- 
gether with ASCI {YMR116C) and HXKl {YFR0S3C), 
whose deletion was shown to induce protracted fermenta- 
tion (Additional files 1 and 4). Hxklp, which is associated 
with heat shock response, is one of three protein kinases 
(Glklp, Hxklp, and/or Hxk2p) required for glucose phos- 
phorylation, as a means of intracellular signalling via 
Gprlp-Gpa2p. Whilst Asclp, a guanine dissociation in- 
hibitor of Gpa2p, links glucose metabolism to nitrogen 
metabolism, as it represses Gcn4p when amino acids are 
sufficient [66]. 



PLCl (YPL268W) was also identified in the FEG data- 
set (Additional files 1 and 4). The inability of the Aplc 
deletant to complete sugar catabolism, as shown in 
Figure 2 alludes to the importance of Plclp in cellular 
metabolism; although the physiological roles during fer- 
mentation are not well understood. Plclp has been im- 
plicated in nutrient induced signal transduction via the 
Gprlp-Gpa2p mediated GPCR system, coupling glucose 
and nitrogen signalling at the sensor level [64,67], as well 
as controlling the MAP kinase cascade via activation of 
Pkclp. This cascade is essential to adaptation to high 
osmolarity and cell wall integrity [41]; obvious important 
attributes for successful fermentation. 

HRKl {YOR267C) and PTK2 {YJR0S9W) were also 
identified and encode protein kinases responsible for the 
activation of Pmalp in response to glucose [68]. The 
membrane potential generated by Pmalp is modulated 
by the uptake of potassium (K^) by the Trklp and Trk2p 
transporters; allowing for different voltage sensitive trans- 
porters to function in cation transport [48,68]. The central 
role of pH (and cation) homeostasis in fermentation is fur- 
ther supported by the inclusion of TRKl {YJL129C) in the 
FEG dataset, which encodes the high affinity K"^ trans- 
porter and is also found to be important for resistance to 
high-glucose [4] and high acidity [25]. 

Interestingly, Pmalp activity is reduced in yeast in 
high sugar fermentation conditions with recovery to nor- 
mal activity being incomplete upon transfer to low sugar 
conditions [4]. Teixeira and co-workers [4] propose that 
vesicle trafficking to the plasma membrane is affected, 
given the mis -targeting of Pmalp in vma mutants lack- 
ing V-ATPase subunits to the vacuole [55]. pH and ion 
homeostatic mechanisms via V-ATPase are implicated in 
both osmotic shock from salt and sugar [4], with V- 
ATPase contributing to not only vesicle protein trafficking 
but intracellular (cytosolic) pH through the translocation 
of protons into the vacuolar lumen. NHXl {YDR456W} 
encoding the endosomal IC (Na"^)/H"^ exchanger, Nhxlp, 
was also identified as a FEG in our study. Nhxlp is impli- 
cated in regulating vacuolar and cytosolic pH [55], and is 
thought to represent an early response to hyperosmotic 
conditions by the transport of ions into the vacuole [69]. 

The complex interaction between Pmalp, V-ATPase 
and other symporters is in part also mediated by the 
lipid rafts associated with these membrane bound pro- 
teins. In a recent study [70] lipid raft integrity of the 
plasma membrane was implicated in pH homeostasis 
and growth. Edelsofine, an anti-tumor lysophosphatidyl 
choline analog, was shown to alter the organisation of 
the plasma membrane, resulting in ubiquitination and in- 
ternalisation of Pmalp to the vacuole, resulting in cyto- 
plasmic acidification [70]. The link between membrane 
(lipid raft) integrity, intracellular pH, and fermenta- 
tion phenotype is also inferred from the sensitivity of 



Walker et al. BMC Genomics 2014, 15:552 
http://www.bionnedcentral.conn/1471 -21 64/1 5/552 



Page 11 of 1 7 



the following FEG mutants to Edelsofine, Aascl, Adfb2, 
Asnf4, AsinS, Atrkl, Avma2, Avma4, AvmaS, Avma6, 
Avma7, AvmalO, Avmall, Avmal6 and Avph2 [70]. 

Role of vacuolar H"^ATPase in fermentation 

As mentioned, 18 genes associated with vacuolar acidifi- 
cation, specifically the vacuolar ATPase complex were 
identified in this study (Additional file 4 and Table 3). 

The vacuolar ATPase complex is the major contribu- 
tor in regulating organelle pH, specifically between the 
cytosol and vacuolar lumen [49] . This action is intrinsic- 
ally connected to cellular ion homeostasis (through the 
action of anti-porters utilizing the generated membrane 
potential via a proton electrochemical gradient), as well 



as proteolytic activation (digestion of autophagasomes 
during autophagy or vacuolar proteolysis induced by 
starvation), endocytosis, and vacuole fusion [49,71]. Fur- 
thermore, intracellular pH and Ca^^ ion concentrations as 
well as direct interaction of the V- ATPase complex with 
the actin cytoskeleton are implicated in bud formation 
and polarized growth during cell cycle [50,72]. In relation 
to the vacuoles lytic function, VAM3 {YOR106W) encod- 
ing the syntaxin-like protein VamSp (a vacuolar t-SNARE) 
is responsible for the correct trafficking and processing of 
proteinases A, B and carboxypeptidase Y and maturation 
of alkaline phosphatase [73]. The relevance of VamSp to 
fermentation outcome and tolerance of stress (hyperos- 
molarity, heat, ethanol and oxidation) is highlighted by the 



Table 3 Effect upon fermentation of yeast lacking genes encoding vacuolar ATPase complex and associated 
proteins 



Gene 


V-ATPase subunit 
Vi domain 


Function (eight subunit peripheral domain) 


Essential for 
fermentation 


VMM {TFPl) 


A 


ATP hydrolysis 


Yes 


VMA2 


B 


Regulatory, ATP binding?, Actin binding? 


Yes 


VMA5 


C 


Stator 


Yes 


VMA8 


D 


Rotor 


Yes 


VIVIAN 


t 


Stator 


Yes 


I /A /! /\ 7 


r 

r 


Rotor 


Yes 


\/l\AA 7 n 




OLdLUl 


T cb 


VMM 3 


H 

Vo domain 


Stator 

(six subunit integral domain) 


Yes 


VP HI 


a (vphl p) 


Proton pore, stator, sorting (vacuole) 


Yes 


STVl 


a (Stvl p) 


Proton pore, stator, sorting (golgi/endosome) 


No 


VMA3 {CUPS) 


c 


Proton pore, rotor (dicyclohexylcarbodiimide binding) 


Yes 


VMAll {TFP3) 


c' 


Proton pore, rotor 


Yes 


VMA16{PPA1) 


c" 


Proton pore, rotor 


Yes 


VMA6 


d 


Rotor 


Yes 


VMA9 


e 

Assembly Factors 


? 


nd 


VPH2 {VMA12) 


Vph2p 


Integral membrane protein; vacuolar ATPase (V-ATPase) assembly 


Yes 


VMA22 


Vma22p 


Peripheral membrane protein; vacuolar ATPase assembly 


Yes 


VPS3 


Vps3p 


CO RVET tethering complex; cytoplasmic protein required for sorting & processing of soluble 
vacuolar proteins, acidification of vacuolar lumen, & assembly of V-ATPase 


No 


PKRl 


Pkrlp 

RAVE complex 


V-ATPase assembly factor, functions with other V-ATPase assembly factors in ER to assemble 
V-ATPase membrane sector (Vq) 


Yes 


SKPl 


Skplp 


Evolutionarily conserved kinetochore protein; part of SCF ubiquitin ligase complex, CBF3 
complex binding centromeric DNA, & RAVE complex regulating assembly of V-ATPase 


nd 


RAVI 


Ravi p 


Subunit of RAVE complex (Ravlp, Rav2p, Skplp), promotes assembly of V-ATPase 


Yes 


RAV2 


Rav2p 


Subunit of RAVE complex, promotes assembly of V-ATPase 


Yes 



Fermentations (100 mL) were conducted in CDGJM_200 as described in Materials and Methods. Deletants which resulted in protracted fermentation were 
considered essential for fermentation in high sugar media, whilst those which were not affected, were considered non-essential. 
VI\/IA9 {YCL005W-A) and SKP^ {YDR328C) are not in the homozygous diploid deletion library. SKP^ is an essential gene for growth, 
nd: not determined. 
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protracted fermentation of the AvamS mutant (Additional 
file 1) and growth sensitivity towards these stress parame- 
ters (Additional file 2). 

The importance of the assembly and function of the 
V-ATPase complex during fermentation, was highlighted 
by the number of genes identified (Additional files 4 and 
Table 3). These encode either the subunits of the Vi per- 
ipheral membrane and Vq integral membrane domains, 
or the assembly factors including the RAVE complex, 
which are required for the reversible assembly and disas- 
sembly into the free Vi and Vq domains, in response to 
glucose [50,51]. The assembly of V-ATPase is thought to 
prevent cytosolic acidification when metabolism is ac- 
tive, as protons are sequestered in the vacuolar lumen 
by the V-ATPase. This is particularly relevant in high 
sugar fermentation, when the plasma membrane bound 
pump Pmalp, the principal regulator of cytosolic pH, is 
transiently inhibited [4]. As glucose becomes limited, the 
disassembly of some V-ATPase complexes is thought to 
conserve ATP, as ATP hydrolysis is absent in the isolated 
Vi domain. The role of vacuolar ATPase in fermentation 
is supported by previous research showing deletion of 
structural genes resulting in growth sensitivity to ethanol 
[8-10] and hyperosmolarity due to high concentrations 
of glucose (300 g L"^, [4]). How these differences relate 
to the pH-dependent growth of vma (vacuolar mem- 
brane ATPase) deletants [50] is unclear. Whilst vma 
deletants, through the loss of V-ATPase activity, exhibit 
growth sensitivity at high pH (>pH 7) or high extracellu- 
lar calcium and zinc, and reportedly are unable to grow 
on non-fermentable carbon sources, they are able to 
grow under acidic conditions (pH 5.5) albeit slower than 
the wild type [50,74]. In this study, the mutants were 
able to grow in CDGJM_200 (pH 3.5) after an extended 
lag period yet only reached 60-80% of the parent bio- 
mass (Figure 5). The ability of these cells to grow even 
in the presence of initially high sugar and acid and later, 
increasing ethanol, is strong evidence that vacuolar 
acidification by this enzyme complex is not the only 
mechanism involved. Yeast possess additional, and yet 
unidentified independent mechanisms to acidify the 
endosomal compartments including the vacuole. We have 
preliminary microscopy data using 6-carboxyfluorescein 
diacetate which supports this notion. Specifically, early in 
fermentation (24 h) vacuolar acidification, as anticipated, 
was reduced in a Avmal mutant compared to the parent 
BY4743. However, after 72 h and 125 h, by which time fer- 
mentation was progressing well, vacuolar acidification was 
restored in the mutant (data not shown). Recently, va- 
cuolar acidification, particularly the proton dependent 
storage of neutral amino acids rather than protein degrad- 
ation in the vacuole, has also been linked to chronological 
aging [38]. The authors proposed that age-induced mito- 
chondrial dysfunction is a result of decreased vacuolar 



acidification in aging mother cells, and the regeneration of 
vacuolar acidification in budding daughter cells, is coinci- 
dental with rejuvenation of lifespan. This network between 
the vacuole and mitochondria, thought to be mediated by 
nutrient-sensing pathways such as PKA, Sch9 and TOR 
[38], is important given that fermentation occurs over an 
extended time-frame, primarily by stationary phase cells 
which are progressively aging. 

Several proposed mechanisms to maintain vacuolar 
pH independent of V-ATPase include passive diffusion 
and dissociation of weak acids, endocytic internalisation 
of acid equivalents or Pmalp itself [74], and ammonium 
ions acting as protonophores after uptake [74]. These 
mechanisms would enable growth under acidic conditions 
through pH equilibration across the plasma membrane 
and vacuolar membrane, maintaining the intracellular pH 
gradient between the cytosol and vacuole, relative to the 
outside. This process is vital to cell adaption to fermenta- 
tion associated stresses, and allowing for complete sugar 
catabolism albeit at a slower rate. Whilst some of these al- 
ternative mechanisms may play a role normally in pH 
regulation, V-ATPase is the principal proton pump pro- 
posed to protect yeast from heat stress [75], alcohol stress 
[8,9], osmotic stress [4], and acid stress [25] associated 
with fermentation, as well as against stresses of air-drying 
[29]. The latter is important in this context since commer- 
cial wine yeast is generally supplied as active dried yeast. 

Whilst there are a myriad of processes and signalling 
cascades related to stress response, it is only recently 
that the inter-relationships are being elucidated. For ex- 
ample V-ATPase and the high-osmolarity glycerol (HOG) 
response pathway are proposed to function in parallel, as 
adaptive mechanisms to prevent salt toxicity under salt in- 
duced stress [76]. Both pathways are crucial for successful 
sugar fermentation, with not only vma mutants being 
identified but also the genes for Hoglp and its activator, 
Pbs2p. It is possible that communication between such 
processes as vacuolar acidification (V-ATPase), protein 
trafficking and turnover, and the TOR and HOG signalling 
pathways is through inositol-containing lipids in stress re- 
sponse signaling [77], in response to stress adaptation. 

Microautophagy and the EGO/GSE complex and 
fermentation 

MEHl and SLM4 encode 2 of the 3 subunits of the EGO 
complex and deletion of either severely affected fermen- 
tation in this study (65% slower than the parent). With 
the third component, Gtr2p, the EGO complex is a va- 
cuolar membrane-associated protein complex required 
for activation of microautophagy [78] through its func- 
tion as an activator of the nitrogen regulatory TORCl 
(Torlp or Tor2p-Koglp-Lst8p-Tco89p) complex [79]. In 
this study TC089 was also identified as a FEG, with 
Atco89 increasing fermentation duration by 43%. Mehlp 
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and Slm4p are also members of the GSE (GAPl Sorting 
in the Endosomes) complex, a GTPase complex required 
for intracellular sorting of Gaplp out of the endosome, 
for eventual delivery to the plasma membrane [80]. The 
GSE complex is also composed of two small GTPases 
(Gtrlp and Gtr2p) and Ltvlp, which do not have a crit- 
ical role in fermentation (Additional file 1). These find- 
ings suggest that MEHl and SLM4 have a major role in 
fermentation, through vacuolar acidification (Mehlp), 
probably through amino acid uptake into the vacuole, 
microautophagy and protein sorting of Gaplp. 

Microautophagy and the EGO/GSE complex have 
been previously implicated in mechanisms allowing yeast 
cells to grow under low temperature and high pressure 
[7]. Abe and Minegishi [7] proposed that the EGO/GSE 
complex may contribute to cell surface delivery of amino 
acid permeases under such conditions, given the marked 
defect in histidine, leucine and lysine uptake in the EGO/ 
GSE complex mutants (Agtrl, Agtr2, Amehl (Aegol) and 
Aslm4 (Aego3)), upon high pressure and low temperature 
incubation. 

As mentioned, the EGO complex is involved in nitro- 
gen sensing via its activation of the TORCl complex, 
such that EGO monitors intracellular levels of leucine 
and glutamine (reviewed in [79]). The inability to acti- 
vate TORCl during fermentation may result in fermen- 
tation arrest via several processes in which TORCl is 
involved: the regulation of amino acid uptake, early gly- 
colysis or induction of autophagy (reviewed in [79]). 

It is known that autophagy is triggered by nutritional 
stress, specifically carbon and nitrogen depletion, and is 
induced as a prelude to autolysis in yeast conducting the 
secondary fermentation of sparkling wine production 
[81]. Autophagy is critical for fermentation outcome in 
primary fermentation as shown in our study, with 13 re- 
lated genes identified as FEGs. In a recent study, Piggott 
et al. [82] looked at competitive fitness (growth) during 
a 12-day fermentation of a deletion library pool (S288c) 
in synthetic grape juice. Autophagy and ubiquitination 
were enriched in the dataset with respect to reduced fit- 
ness. However, our datasets for these GO terms are very 
different, with only 3 genes in common for autophagy; 
{ATG7 (YHRiyiW), SNX4/ATG24 {YJL036W), CIS/ 
ATG31 (YDR022C)). Other contradictions were the obser- 
vations for Adoa4 and Apexl. Deletion oiDOA4 (homozy- 
gous diploid) did not result in reduced competitive fitness 
[82], contrary to our findings, wherebyAdoa4 resulted in 
stuck fermentation. Another anomaly was the finding that 
Apexl (protracted fermentation in our study) resulted in 
increased fitness [82], implying that fermentation would 
be either similar to the parent or faster. The recent revela- 
tion that autophagy occurs not only during starvation but 
very early in fermentation [82], in response to as yet un- 
known stresses, alludes to the importance of the various 



autophagic responses and would explain the observed 
protracted fermentation phenotype in such mutants in 
our study. 

It is evident that a number of mechanisms involving a 
myriad of genes are involved in the adaptive response to 
the major stressors encountered during high sugar fer- 
mentation as such only some of these have been dis- 
cussed to in this study. 

Conclusion 

This study reports on the identification of 93 genes 
whose deletion leads to fermentations that are pro- 
tracted (84 genes) or arrested ("stuck"; 9 genes). To- 
gether, the 93 genes are referred to as Fermentation 
Essential Genes, and are representative of the Termen- 
tome'; a dataset of genes (from laboratory yeast) which 
modulate fermentation. The importance of the PEG is 
further supported through their presence in two fermen- 
tation relevant databases, LIT_FRYG and SGD_FRYG 
(Additional file 3), collated as part of this study. These 
datasets are by no means complete, but will be expanded 
upon as new datasets become available. 

In this study we have built upon the understanding of 
a group of key genes within fermentation. We have iden- 
tified several cellular processes essential for response 
and adaptation to a physiologically stressful environment 
in yeast, namely maintenance of vacuolar acidification, 
microautophagy and sugar related signalling. Whilst in- 
dividual biological processes are alluded to with respect 
to fermentation, the question remains; how do these dif- 
ferent sensors, signals and cellular responses interact 
with each other? A systems biology approach allows us 
to look at what appear to be unrelated processes, to- 
gether as a 'network' or 'matrix', in a holistic approach, 
as well as the traditional reductionist approach [2,19]. By 
taking such a perspective on wine fermentation, we have 
gained a better understanding of a fundamental bio- 
chemical pathway. This understanding allows design of 
new industrial yeast strains better suited to specific in- 
dustrial processes. For example demonstration for the 
importance of ion homeostasis and intracellular pH 
maintenance raises the possibility that Directed Evolu- 
tion could be used to generate new strains using a low 
pH environment as the selective pressure. 

Methods 

Yeast strains and media 

This study used the collection of yeast diploid deletion 
strains developed by the Yeast Genome Deletion Project 
[12]. The homozygous {Aorf::KanMX4) gene deletions 
were harboured in BY4743 {MATa/a, hisSM/hisSM, 
leu2A0/leu2A0, lys2A0/LYS2, METlS/metlSAO, uraSAO/ 
uraSAO (BY4741/BY4742)). Strains were maintained on 
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YPD (1% yeast extract, 2% bactopeptone, 2% D-glucose) 
containing 200 mg L"^ geneticin (G418 sulfate; Amresco). 

Chemically defined grape juice medium containing 
200 g L"^ sugar as equimolar amounts of glucose and 
fructose (high sugar; CDGJM_200) and 450 mg L"^ FAN 
as a mixture of amino acids and ammonium chloride 
[17,45] was used. The CDGJM_200 was supplemented 
with 150 mg L"^ uracil, and included 3 g L"^ polyphenol 
extract (Cat: Tppr, OenoProd, Sari), which was dissolved 
in 100% ethanol (3 g powder in 5 mL, 24 h, in dark) prior 
to being added to the sterile medium. The CDGJM _200 
was stored in the dark for a maximum of 24 h (4°C) before 
use. CDGJM_20 was identical to CDGJM_200 except that 
the sugar concentration was decreased to 20 g L"^ (low 
sugar; CDGJM_20). 

Genome-wide screening for yeast deletants with 
protracted fermentation in high sugar media 
(CDGJM_200) 

Liquid handling of the samples was performed using a 
CAS3800 robot (Corbett Robotics, Sydney). The library 
was replicated (6 (iL) into 0.6 mL YPD in 96 deep-well 
plates (Cat: P-DW-20-C, Pacific Lab Products) and incu- 
bated 24-48 h at 28°C in a humidified box without shak- 
ing. Each plate had 11 rows (Al-lltoHl-11) containing 
single samples of the deletion strains and four replicates 
each of the parental strain, BY4743 (A12-D12) and un- 
inoculated controls, CDGJM_200 (E12-H12). Cells were 
resuspended by agitation on a mini-vortexing-shaker 
and 6 [iL was inoculated into 0.6 mL CDGJM_200 in 96 
deep-well plates and covered with breathable sealing film 
(Cat: BF-400, Adelab Scientific). The initial fermentation 
screen was performed in duplicate in humidified boxes. 
After 180 h at 28°C, the cells were pelleted by centrifuga- 
tion (1620 X g, 5 min) and 0.2 mL supernatant samples 
were dispensed into 96 well flat bottomed plates (Cat: P- 
96-300 F-PS, Adelab Scientific) and frozen at -20°C for 
subsequent sugar analysis. Residual glucose and fructose 
were quantified enzymatically using a commercially avail- 
able kit (Cat: 139106, Boehringer Mannheim) adapted for 
96-well plates. The total volume of the reaction was re- 
duced from 1.62 mL to 0.125 mL. Three times the con- 
centration of the reaction buffer was used to provide 
NADP and ATP to excess, and to ensure accurate pipet- 
ting, 10 (iL of a 1:10 dilution of the enzymes in MQ water 
was used (as opposed to 1 (iL). The samples were diluted 
1:10 (for < 20 g L"^) and 1:100 (for > 20 g L"^) and deter- 
mined against a set of identically diluted standards which 
were within the linear range of the reaction (0-0.08 g L 
~^). Absorbance at 340 nm was measured using a (iQuant 
microplate spectrophotometer (Bio-Tek Instruments). 

A second experiment ((/initial screen) was conducted, 
using a sacrificial plate method [83] on the 336 deletion 
strains identified as having protracted fermentation. 



At specific time points during fermentation (41, 89, 
113, 137 and 190 h), a replicate plate was sacrificed 
to measure growth (ODeoo) and residual sugar in the 
CDGJM_200 media. Absorbance was measured using 
a microplate spectrophotometer. 

Low sugar fermentations were also conducted on the 
336 deletants (cf second screen) using CDGJM_20, in 
which the sugar content was reduced to 20 g L"^ and 
the duration of the fermentation was 56 h. Growth and 
fermentation progress was monitored as previously de- 
scribed at two time points: 32 and 56 h. 

Evaluation of identified yeast deletants in laboratory scale 
(100 mL) fermentations in CDGJM_200 

Triplicate 100 mL fermentations were performed in 
CDGJM_200 using 250 mL Erlenmeyer flasks fitted with 
airlocks [84] . N2 sparging was omitted as previous studies 
demonstrated that anaerobic conditions were achieved 
within 4-6 hours of fermentation commencement (data 
not shown). Fermentation progress was monitored by re- 
fractive index (°Brix). Fermentations were considered dry 
and terminated when the residual sugar was less than 
2.5 g L"^ as determined by Clinitest™ (Cat: 2107, Bayer). 
Residual glucose and fructose were determined by enzym- 
atic analysis as described above. 

Confirmation of the identity of yeast deletants 

Gene deletions in the 93 FEG (Fermentation Essential 
Genes) yeast clones were confirmed by PCR amplifica- 
tion of deletion cassettes from isolated genomic DNA 
using the protocol described in the Saccharomyces Gen- 
ome Deletion Project (http://www-sequence.stanford. 
edu/group/yeast_deletion_project/deletions3.html) . The 
PCR products of select clones were verified by DNA se- 
quencing (AGRF, Adelaide). 

Classification and numeric enrichment of identified genes 
annotated to Gene Ontology (GO) terms using 
computational software tools 

The datasets were analysed using GOToolBox; specifi- 
cally GO-Stats (http://genome.crg.es/GOToolBox/) [30] 
and GO Finder Version 0.83 software (http://www.yeast 
genome.org/cgi-bin/GO/goTermFinder.pl) which allowed 
the hierarchical clustering and over-representation (en- 
richment) or under-representation (depletion) of genes 
based on shared Gene Ontology (GO) terms. Both com- 
putational tools use a hypergeometric distribution with 
Multiple Hypothesis (Bonferroni) Correction (and add- 
itional False Discovery i^ate; GO Finder) to calculate p- 
values. The clustering of genes based on their shared gene 
annotations was also determined using GO-Proxy (http:// 
genome.crg.es/GOToolBox/) and/or SGD Gene Ontology 
Slim Mapper (http://www.yeastgenome.org/cgi-bin/GO/ 
goSlimMapper.pl) . 
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The interpretation of the GO terms identified was 
simpUfied using the web-based visuaUsation tool GO- 
Module (http://lussierlab.org/GO-module) [31]. GO 
biomodules are 'nodes' which represent the key GO 
terms (K nodes) and their hierarchical descendents which 
are considered true-positive GO terms (T nodes). False 
positive prioritized GO terms in the input data are identi- 
fied as F nodes. 

Two databases denoted Fermentation Relevant Yeast 
Genes (FRYG), were compiled from data available on 
the Saccharomyces Genome Database (http://www.yeast 
genome.org/) and from a literature survey. The data- 
bases referred to as SGD_FRYG and LIT_FRYG were 
used in comparison studies. Statistical analysis of the re- 
sults was performed using a hyper geometric algorithm 
as described by Martin et al. [30] to determine whether 
there was an enrichment (over representation) of genes 
identified in this study when compared to the FRYG 
datasets. 

Availability of supporting data 

The data set supporting the results of this article is avail- 
able in the LabArchives repository, under Fermentation 
Microbiology doi:10.6070/H4057CWl and hyperlink to 
dataset(s) in https://mynotebook.labarchives.com/share/ 
Fermentation%2520Microbiology/MjQuN3wOMDY2Ni 
8xOS9UcmVlTm9kZS8xMzY50DkyMTYxfDYyLjc=. 

The data set supporting the results of this article is in- 
cluded within the article (and its additional file(s)). 

Additional files 



identified. Deletion of these genes resulted in protracted fermentation 
(20% and greater when compared to parent, BY4743). See Additional 
file 1 for fermentation duration values. Two complementary databases of 
yeast genes related to fermentation, referred to as FRYG (acronym for 
Fermentation Relevant Yeast Genes) were generated. The first (LIT_FRYG) 
was compiled from relevant literature, examining single stress conditions 
relevant to fermentation. The second (SGD_FRYG), a compilation of 
phenotype terms as accessed on the SGD, whereby only null mutants 
(deletants) with a specific phenotype response, i.e. decreased resistance 
were selected. The 93 FEG were compared to the two FRYG datasets. 
Genes unique to this fermentation study (absent in the LIT_FRYG and 
SGD_FRYG) and two specific subsets of genes within the 93 FEG dataset 
("Stuck mutants" and "vacuolar (H^) ATPase") are shown. 

Additional file 4: Classification and numerical enrichment of the 
FEG annotated to biological process, function, and cellular 
component. Genes were classified using the computational tools Slim 
Mapper (SGD database; (http://www.yeastgenome.org/) and GOToolBox 
into categories based on GO terms for biological process, function and 
cellular component. Descriptions of the 93 FEG are taken from SGD 
database. Enrichment of genes annotated to particular GO terms includes 
corresponding p values. GO terms were organised into GO biomodules 
using the web based visualisation tool GO-Module whereby K nodes refer 
to the key GO terms and T nodes refer to hierarchical descendents of K, 
regardless of statistical value (p value). F nodes represent the false 
positive prioritised GO terms in the input data. F nodes determined for 
particular GO terms were not precluded because of the nature of the 
GO term. Genes annotated to GO terms not defined through GO module 
and not significantly enriched were classified into groups using SGD Slim 
Mapper. 
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